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O ■ Abstract 

$h ■ We consider cosmological models with a scalar field with equation of state w > 1 that contract 

towards a big crunch singularity, as in recent cyclic and ekpyrotic scenarios. We show that chaotic 
mixmaster oscillations due to anisotropy and curvature are suppressed, and the contraction is 



X 



described by a homogeneous and isotropic Friedmann equation if w > 1. We generalize the results 
to theories where the scalar field couples to p-forms and show that there exists a finite value of w, 



(N 
> 
On 

depending on the p-forms, such that chaotic oscillations are suppressed. We show that Z2 orbifold 

(N 

compactification also contributes to suppressing chaotic behavior. In particular, chaos is avoided 

m ; 

. in contracting heterotic M-theory models if w > 1 at the crunch. 
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I. INTRODUCTION 



In cosmological models with a big crunch/big bang transition, a key issue is the behavior 
as the universe contracts towards the crunch. From the classic studies of Belinskii, Khalat- 
nikov and Lifshitz (BKL) |2J, |3( and others , it is known that the contraction 

can either proceed smoothly or chaotically. These studies have focused on models in which 
the universe contains matter and radiation, or, more generally, an energy component whose 
equation of state is w < 1 (where w = p/p is defined as the ratio of the pressure p to the 
energy density p). If w < 1, a contracting homogeneous and isotropic solution is unstable to 
small perturbations in the anisotropy and spatial curvature. As the overall volume shrinks, 
the anisotropy causes the universe to expand along one axis and contract along the others, 
a state that can be approximated by the anisotropic Kasner solution. The spatial curvature 
causes the axes and rates of contraction to undergo sudden jumps from one Kasner-like 
solution to another, an effect known as "mixmaster" 0, n| behavior. If the curvature is 
not spatially uniform, then the chaotic behavior in different regions is not synchronized and 
the universe becomes highly inhomogeneous at the big crunch. Hence, mixmaster behavior 
could potentially wreak havoc in cosmological models with a big crunch/big bang transition, 
making them inconsistent with the observed large scale homogeneity of the universe. 

In this paper, we show that the behavior of the universe as it approaches the big crunch is 
very different if there is an energy component with w > 1. The chaotic behavior is suppressed 
and the universe contracts homogeneously and isotropically as it approaches the singularity. 
The reason is that the anisotropy and curvature terms in the Einstein equations grow rapidly 
and become dominant if w < 1, but they remain negligible compared to the energy density 
if w > 1. In the latter case, the Einstein equations converge to the Friedmann equations 
with purely time- dependent terms, a condition sometimes referred to as "ultralocality." The 
effect can be viewed as a generalization of the "cosmic no-hair theorem" invoked in a rapidly 
inflating universe. Here we demonstrate analogous behavior in a slowly contracting universe 
with w > 1. A related result of Dunsby et al [ll, 12, 3] shows that models with < w < 1 
but with p 2 terms in the stress-energy tensor are also driven towards isotropy. 

The cosmic no hair theorem for a contracting universe containing a perfect fluid with 
w > 1 is discussed in section |H] A common example of a perfect fluid is a scalar field <p with 
a potential V(<f)). In section ITTT{ we consider the interaction of the scalar (f> with a p-form 



2 



field F p+ i through an exponential coupling, 



e r p+li 



(1) 



where A is a constant [14]. We consider this case because scalar fields with exponential 
couplings to p-form fields are common in Kaluza-Klein, supergravity and superstring models. 
For the case w = 1, it is known that the contraction is not chaotic if A lies within a 

bounded interval. Here we show that, for any A and p, there is a critical value w cr i t (X,p) 
for which the chaotic behavior is suppressed if w > w CT i t (X,p). 

Our results are of particular importance for the recent ekpyrotic [15J] and cyclic [16J 
cosmo logical models, which have a big crunch/big bang transition with a contraction phase 



dominated by a scalar field with w > 1 



17| . The evolution of perturbations leading up to 



and passin g th rough the transition is an important aspect that remains unsettled 
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251 ] . and may depend on the precise physical conditions leading up to the 
bounce [26, 27]. The present work may be relevant since it suggests that the universe can 
remain homogeneous and isotropic on large scales. Once the evolution becomes ultralocal, 
the whole universe is following the same homogeneous and isotropic evolution all the way 
to the big crunch. 

In section ITVl we explore how time- variation of w affects our conclusions, and in particular 
how w approaching w cr it from above may suppress chaotic behavior. In section we discuss 
some specific models. In particular, we show how orbifolding can remove p-forms that might 
induce chaotic behavior and discuss the special case of heterotic M-theory, which, to leading 
order in the eleven dimensional gravitational coupling k, is on boundary between chaotic 
and smooth behavior. 



II. A "COSMIC NO HAIR THEOREM" FOR CONTRACTING UNIVERSES 

The cornerstone of the inflationary paradigm is an argument known as the "cosmic no-hair 
theorem" , according to which a universe containing a perfect fluid component with w < — 1/3 
will rapidly approach flatness, homogeneity and isotropy at late times, for a wide range of 
initial data (na mely those for which the space curvature, inhomogeneity and anisotropy are 
not very large) [18]. In the Friedmann equation, the energy density for a component with 
equation of state w is proportional to l/a x , where the exponent x = 3(1 The anisotropy 
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term is proportional to a -6 and the spatial curvature term is proportional to a~ 2 . As the 
universe expands, the contribution with the smallest values of x redshifts away more slowly 
than components with larger values of x and so come to dominate the Friedmann equation 
and the components with the smallest value of x overall ultimately dominate. If the energy 
component with the smallest value of w has w < —1/3, then x < 2 and this component 
dominates. For a wide range of initial data, convergence to a homogeneous and isotropic 
expanding universe is assured. 

Below, we will present an analogous "cosmic no-hair theorem" for contracting universes. 
In a contracting universe, the component with the largest value of x will dominate the 
Friedmann equation. Starting from an inhomogeneous and anisotropic initial state, we will 
show that the existence of a perfect fluid with w > 1 (or x > 6) will suppress chaotic 
behavior, and enable a smooth and isotropic contraction to the big crunch. We will find 
that curvature plays a more complicated role compared to the case of expansion. Hence, 
we first obtain a cosmic no-hair theorem for the case of zero spatial curvature and then 
generalize to the case of arbitrary spatial curvature. We intentionally take a pedagogical 
approach that encompasses known results for w < 1 to make our discussion self-contained. 
Our analysis assumes the initial inhomogeneity is small; it is possible that the universe 
evolves towards other attractors for sufficiently large deviations from homogeneity. Our 
conventions are given in [l^ . 

All of our computations are performed in synchronous gauge, 

ds 2 = -dt 2 + h ab (t, x) dx a dx\ (2) 
where we use our freedom to choose a spatial slicing to ensure that the big crunch occurs 



everywhere at t — (det h a b — > as t — > 



state p = wp, the Einstein equations are 



.i 



For a perfect, co moving fluid with equation of 



— k/ + Kj ^ = - — Jp, (3a) 
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where P a is the Ricci tensor on spacelike surfaces, and K a t> is defined by 

1 9 u I A \ 

Kah = 2dt 

= K aj h jb . (4b) 

Near the big crunch, the dynamics of the metric (j2J are ultralocal j^l Q, 0, 13 ■ That is, 
the evolution of adjacent spatial points decouples because spatial gradients increase more 
slowly than other terms in the equations of motion. Therefore, analyzing the dynamics of 
this metric near the singularity and at fixed spatial coordinate x is equivalent to analyzing 
the much simpler system 

ds 2 = _ d t 2 + J2 e 2 ^ x °V%; xo)a<%; x ), (5) 

ij 

where the are y-dependent one-forms that are linearly independent at each point and 
form a homogeneous (but possibly curved) space such as Bianchi type IX [l(3|. The 
which do not depend on y, describe the (generally anisotropic) contraction of this space. 
Both the a^ 1 ' and the (3ij depend on the parameter xo, the spatial point being studied. The 
dynamics of the inhomogeneous universe at a fixed spatial point can be approximated, near 
t = 0, by the dynamics of a homogeneous (but curved and anisotropic) universe. Differences 
in curvature and anisotropy between different x are encoded in the different and Pij 
associated with these points. 

In each Kasner-like epoch, we may perform a rotation so that (3 is diagonal. Furthermore, 
we may separate out the trace of (3 and write it as the "volume scale-factor" a(t), in analogy 
to the isotropic Friedman-Robertson- Walker universe, to obtain the metric 

ds 2 = -dt 2 + a 2 (t) e mt) (a^) 2 , (6a) 

i 

/3i(t)+(3 2 (t)+/3 3 (t) = 0, (6b) 

where the dependence of a(t), the Pi and the on xo has been suppressed. The com- 
bination oe" can be thought of as the effective scale factor along the i th direction, and 
the functions /3j then describe the contraction or expansion of each direction relative to the 
overall volume contraction. We may use our freedom to rescale the a to ensure that at some 
time t , a(to) = 1, A (to) — and det(a^\ a^ 2 \ a^) = 1. Quantities with a subscript zero 
(such as p ) refer to their values at this fixed time. 



The Einstein equations Q close with the equation of energy conservation for the fluid, 

^ = -3(l + u>). (7) 
a log a 

For constant w, this equation has the familiar solution, 

p(a) = p a-^ +w \ (8) 

While we could have included several perfect fluids, with different equations of state Wi, the 
fluid with the largest equation of state will always dominate near the crunch, so it is sufficient 
to consider only one energy component. We have taken this fluid to be co moving, because 
small perturbations of a comoving background are suppressed in a w > 0, contracting 
universe. In particular, the T°j terms that would appear on the right hand side of (j3b|) grow 



only as t 2 /( 1+w ) ; which is slower than the t 2 rate at which the diagonal terms grow 29 1. 



A. The Curvature Free Case 

We first examine the case of Ricci flat spatial 3-surfaces, for which P a b = 0. In this case, 
we write erW = dx\ Then, the Einstein equations reduce to 

^ -\{Pl + Pl + Pl)=P, (9a) 

Pi + 3-A = 0, (9b) 
a 



where a dot indicates a derivative with respect to the proper time t. Integration of (J9bj) 
gives, 

Pi = Qa" 3 , (10) 

while the constraint ()6bjl implies, 

ci + c 2 + c 3 = 0. (11) 



Combining these results, equation (j9a|) becomes a Friedmann equation, 

3 - = pia) + — - = + — , (12) 

\aJ a 6 a 3 ( 1+w ) a 6 

where we define 

° 2 = \{c\ + cl + cl). (13) 
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An anisotropic universe has j3i 7^ 0, i.e. q 7^ 0. The constant a 2 parameterizes the 
anisotropic contribution to the Friedmann equation in (|12|). The anisotropy evolves as 
1/a 6 or x = 6. We define the fractional energy densities fl p and Q a as 

^ = i VT^l e » ( 14a ) 

o 2 /a 6 

^ = ( rr 2/ 6 - ( 14b ) 

These quantities represent the contribution of the perfect fluid and anisotropy to the critical 
density for closure of the universe. Since we are neglecting curvature, fl p + Vt a = 1. 
The solution for the /3j as a function of the scale factor a is, 

P l (a) = c t V3 f ^{pW + a 2 )' 112 . (15) 

J a a 

The limits of integration have been chosen to ensure /3j(l) = 0. For the remainder of the 
paper, we will assume a universe contracting towards a — > as t approaches zero from below. 
Let us now examine the behavior of these solutions for various w. 



1. w < 1: 

When w < 1, the p(a) part of the integral (fTB]) is negligible as a — > 0, and so the solution 
converges to the vacuum (p = 0) Kasner universe during contraction, 

/ t \ 1/3 

a(t) = - (16a) 



to 

A(t) = ^=lnf^). (1Gb) 



cr\/3 ^0 

The Kasner universe is parameterized by three Kasner exponents pi, 

p - 4 + ^ < 17) 

The scale factors in (jBJ) are powers of t: 

ae ft = \t/t \*, (18) 

and the relations (fTTj) and (jlHj) become 

Pi +P2 +P3 = 1 (19a) 
Px+pl+Pl = h (19b) 
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known as the Kasner conditions. These describe the intersection of a plane, the Kasner 
plane, and a unit sphere, the Kasner sphere, as illustrated in Fig. [TJ We will denote the 
intersection, which represents the allowed values of the Pi, as the Kasner circle. The out- 
ermost circle in Fig. corresponds to the limit where w < 1, as the energy density scales 
away and only a vacuum, anisotropic universe remains. 

There are three degenerate solutions where exactly one of the pi is one, and the other 
exponents are zero (the solid black circles in Fig. Q). At all other points on the (dashed) 
Kasner circle exactly one of the Pi is negative. Thus, although the geometric mean of the 
three scale factors a(t) = \t\ is contracting, a single scale factor corresponding to the negative 
Kasner exponent is undergoing expansion to infinity. 

For the curvature-free case, the universe becomes increasingly anisotropic near the big 
crunch if w < 1. In particular, the isotropic solution, p\ = P2 = P3 = 1/3, is inconsistent 
with the Kasner conditions ()19ap and ()19b|) . 



Inspection of (JBJ) reveals that, when w = 1, the matter density and the anisotropy terms 
in the Friedmann equation (fT2"j) scale with the same power of a, so Q p and remain fixed. 
The solutions are 



This solution is very similar to the p = case, and indeed we may define the Kasner 
exponents, 



2. w = I: 




(20b) 



(20a) 




(21) 



The Kasner conditions are different. If we define 




3 a 2 + p 3 



(22) 



then the Kasner conditions are 



Pi +P2 +P3 = 1, 



(23a) 



pl + v\ + pi = i - q 2 




(23b) 
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FIG. 1: The Kasner plane pi + P2 + P3 = 1 and its intersections (the Kasner circles) with various 
spheres p\ -\-p\ + p 2 = 1 — q 2 where q 2 = |(1 — f2 CT ); see (|22|l. The vacuum solution corresponds to 
£l a = 1 (the outermost circle). The inner circles are relevant to the case where w = 1 and f2 CT < 1. 
In the white regions, the Kasner exponents are all positive (corresponding to contraction); in gray 
regions, one exponent is negative (expanding). If the spatial curvature is non-zero, points along the 
circles in the white region (thick parts of circles) are stable but points in the gray regions (dashed 
parts of circles) are unstable, jumping to new values after a short period of contraction. If a model 
(i.e. a circle) has an open set of stable points (the three innermost circles but not the outermost 
circle), the contracting phase does not exhibit chaotic mixmaster behavior. 

The first condition is unchanged from (jl9b|) but the right hand side of the second condition 
has been modified. Increasing Q a corresponds to increasing the radius of the Kasner sphere. 

The w — 1 model allows us to explore the behavior of the contracting universe as a 
function of fl a . The perfectly isotropic case corresponds to D, a = 0, which is the usual flat 
Friedmann-Robertson- Walker solution (innermost circle, in the limit where the circle has 
shrunk to a point, in Fig. Q). Unlike the vacuum Kasner case, all of the Kasner exponents 
are positive (i.e. lie within the white region of Fig. ^) provided that Q a < 1/4 (within 
the larger, solid circle inscribed in the triangle). For this range, none of the scale factors 
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is increasing during the contraction, although they are decreasing at different rates. When 
Q a > 1/4 (third largest circle), then some points on the Kasner circle have a negative Kasner 
exponent (dashed part of circle) and other points may have all positive Kasner exponents 
(solid, thick parts of circle). 

Thus, ignoring the curvature, the w = 1 case with non-zero Q a contracts smoothly but 
anisotropically to the crunch. In the special case where Q a = 0, the contraction is isotropic. 



For w > 1, the energy density dominates (Q p — ► 1) as a — > 0, and the metric approaches 
the approximate form 



where we have chosen the constants of integration so f3{ = at t = to. The crucial feature is 
that the time- varying part of the $ is proportional to t a where a is positive if w > 1. This 
means that the approach a constant and the universe becomes isotropic at the crunch 
Q. 

This simple result is a "no-hair theorem" for universes without spatial curvature: When 
w > 1, an initially anisotropic universe becomes isotropic (Q a — ► 0) near the big crunch. 
The w > 1 case is stable under anisotropic perturbations. For w < 1, the universe becomes 
increasingly anisotropic in the sense that Q a —>■ 1 as a —>■ 0. For w — 1, CT remains fixed as 
a — > 0. Evolution is smooth (no mixmaster behavior) in all cases, and is well-approximated 
as a Kasner metric with constant coefficients for sufficiently small a. 

B. Curvature and Chaos 

Complex behavior can arise when there is non-zero spatial curvature in a contracting 
universe. This may seem surprising at first, since the spatial curvature for a homogeneous 
and isotropic universe grows as 1/a 2 , which increases more slowly than either the anisotropy 
or the energy density of a component with w > —1/3. However, we have seen above that 
the contracting phase for w < 1 is anisotropic. We will show below that this can produce 



3. w> 1: 




(24b) 



(24a) 
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rapidly growing curvature perturbations and chaotic behavior. On the other hand, we will 
see that chaotic behavior is suppressed if w > 1 and the contraction approaches isotropy as 
a -> 0. 

We now allow the crW to have an x-dependence and consider a curved manifold. The 
spatial Ricci tensor for the metric (JSJ) has the form [| 



The functions S a b ijk depend only on the a % and their space derivatives, and are independent 
of time. 

The expression (J25|) reveals a crucial connection between the behavior of anisotropy and 
curvature near the big crunch. In the isotropic limit, $ = and (J2"5|) reduces to the 
homogeneous and isotropic 1/a 2 scaling discussed above. However, the terms in (|25|) are 
essentially ratios of scale factors. Thus, if the anisotropy is growing as a — > 0, some terms - 
involving ratios of expanding and contracting scale factors - will grow, and the corresponding 
curvature components will scale faster than 1/a 2 . For w < 1 the anisotropy dominates near 
the crunch, and, as we will discuss below, this causes the curvature to grow and induce 
chaos. By contrast, in the w > 1 model, the anisotropy vanishes at the crunch, and the 
curvature scales as the usual 1/a 2 , which may be neglected. 



In this case, we begin by assuming that the behavior near the crunch is described by 
the vacuum Kasner solution, with Kasner conditions (|19aj) and f!19b|) . Using the Kasner 
solution, it is readily seen that the Einstein equation ()3a|) contains a leading order term 
with time dependence t~ 2 . 

The second Einstein equation (j3b|) is a consistency check for our assumption of ultralo- 
cality. For an appropriate choice of the a® - a basis for one of the Bianchi universes - this 
equation vanishes identically and the metric (jUJ) solves the Einstein equations. 

The third Einstein equation (J3*c|) indicates that the simple Kasner solutions must break 
down near the big crunch. If we order the Kasner exponents as p\ < P2 < P3, then the most 
divergent term in the third Einstein equation comes from terms in the spatial curvature 




(25) 



1. w < 1: 
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(J23J), with leading time dependence, 

t -2(i-2 Pl ) i (26) 



The leading term is more divergent than t~ 2 , since the Kasner conditions (J19a[) and (jl9bj) 
imply pi is always negative. Therefore, our smoothly contracting solutions are not stable to 
perturbations in the spatial 3-curvature. A small amount of curvature will grow and come 
to dominate the dynamics before the big crunch. 

The behavior of the universe in this regime has been extensively studied and is known to 
be chaotic . The spatial curvature terms cause the Kasner exponents pi and 

the principal directions o % to become time- dependent during contraction. 

More precisely, the exponents and principal directions are nearly constant for stretches of 
Kasner-like contraction, during which the curvature is negligible. These Kasner-like epochs 
are punctuated by short intervals when the curvature momentarily dominates. The ex- 
ponents and principal directions suddenly jump to new values, and then a new stretch of 
Kasner-like contraction begins during which the curvature terms are again negligible. The 
universe undergoes an infinite number of such jumps before the big crunch. The chaotic, 
non-integrable evolution is equivalent to that of a billiard ball 8], which experiences free 
motion interrupted by collisions with walls. Models with this oscillatory behavior are called 
chaotic. 

This presents a problem for cosmological models, as one expects curvature perturbations 
in any realistic universe will cause the local value of the curvature to vary from point to point. 
If each spatial point evolves independently and chaotically, the evolution of nearby points 
diverges very quickly as contraction continues, and the universe rapidly becomes highly 
inhomogeneous as a — > 0. If w < 1 throughout the contracting phase, it seems unlikely 
that the observed homogeneous universe could emerge from this state after the bounce to 
an expanding phase [31 1. 

2. w = I: 

The chaotic behavior is mitigated in the w — 1 case. Recalling our discussion of the 
curvature-free scenario, it is clear that there are regions of non-zero measure on the Kasner 
circle for which all of the Pi are positive. We will refer to these points as stable. All choices 
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of pi when Q a < 1/4 are stable. If the universe begins at a stable point, the curvature term 
remains negligible as a —>■ and the contraction is smoothly Kasner-like. 

However, when Q a > 1/4, some choices of the Pi will have one Pi < 0. If the universe 
begins at one of these points, the curvature term will grow and become dominant, causing 
the values of Pi and the principal axes cr, to change. We refer to these points as unstable. A 
more complete analysis |3j reveals that, after a finite number of jumps, the universe hits a 
point in the open set of stable pi. From this point onwards, the universe contracts smoothly 
and without any further jumps. 

We call these models non-chaotic, since the universe is guaranteed to arrive at a stable 
point as a — > 0. Non-chaotic models (Kasner circles) may contain both stable and unstable 
points, but they will always oscillate only a finite number of times before arriving in the set 
of stable points, after which the behavior is integrable. 

3. w> 1: 

For w > 1, curvature does not affect the contraction. The key is the time- dependence 
of the (3i in (|24|) . which approach zero as a positive power of t as t — ► 0. Consequently, the 
exponential factors e*' in the metric approach constants. The leading order time-behavior 
of Pa is simply that of a homogeneous and isotropic universe, 

1 4_ 

p ~ — ~ \t\ 3d+»). (27) 
a 2 

This is always less divergent that t~ 2 for w > 1. Thus, even in the presence of initial 
anisotropy and curvature, the solution for w > 1 converges to the isotropic solution repre- 
sented by the central point on the Kasner sphere in Fig. 1. 

We can generalize our cosmic no-hair theorem (described at the end of section III Aj) to 
include models with spatial curvature. The Einstein equations for a contracting universe 
with anisotropy and inhomogeneous spatial curvature converge to the Friedmann equation 
for a homogeneous, flat and isotropic universe if it contains energy with w > 1, and that 
for a homogeneous, flat but anisotropic universe if w — 1. The w < 1 case becomes highly 
inhomogeneous and the no-hair theorem is inapplicable. 
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III. COUPLING TO P-FORMS AND CHAOS 



In section |Hj we assumed that the evolution of the universe was dominated by an energy 
component with fixed equation of state evolving independently of other matter in the uni- 
verse. The component could have been a scalar field or a perfect fluid. We found chaotic 
behavior for w < 1 in the presence of curvature but non-chaotic behavior for w > 1. 

In this section, we want to consider how the behavior for w > 1 can change if the fluid 
is imperfect or couples to other components. In many theories, including Kaluza-Klein, 
supergravity and superstring models, the relevant energy consists of a scalar field that is 



coupled to p- 



brms. Consequently, we will focus on this important example, as others have 



in the past 



mm 



To determine the effect of the coupling to p-forms on chaotic behavior, our approach is 
similar to our analysis for spatial curvature, where we assume an initial state in which the 
spatial curvature is negligible and then check that it remains small. Here we assume that 
the p-form field strength is initially negligible and ask how its contribution evolves relative 
to the energy density with equation of state w. Our action is 

S = J #xJ=g(\R - \{d<pf - V&) - ^rye^F n+1 2 ), (28) 

where g is the metric, R is the scalar curvature, V is a potential for the scalar field 0, p 
is the rank of the p-form, F is the associated field strength tensor and A is the coupling 
constant. The potential V(<p) is chosen to give fixed equation of state w > 1 in the absence 
a p-form coupling: 



V{<p) = -Voe-V^)*, (29) 

where Vq is a positive constant. Throughout this paper, we assume without loss of generality 
that — > — oo as a — > 0. 

For a given equation of state w and p-form rank, the behavior of the system as t — > 
depends on the coupling A. We can extend the terminology introduced earlier to describe 
the properties for a given A. We classify the p-form coupling parameter A as supercritical 
if the p-form terms grow relative to the scalar field energy density. We call these models 
supercritical, as opposed to chaotic, because if w > 1 it is not known whether chaos occurs or 
whether the p-forms merely play a non-negligible role in integrable dynamics. In the special 
case w = 1, chaos is known to occur, and we call these models chaotic IsL lo, LZl ISll - Values of 
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A for which the contracting solution with negligible p-forms is stable are called non- chaotic 
(some authors use subcritical) . These two cases are analogous to those introduced in section 
im If A is on the boundary between supercritical and non-chaotic, we call A critical. The 
behavior of critical models may be novel, and will be discussed at the end of this section. 

We are assuming that initially the spatial curvature, the anisotropy and the p-form terms 
are small, and then we check if these conditions are maintained as the universe contracts. 
Since we are considering models where w > 1, the model is non-chaotic if the p-forms 
are negligible. The universe may be approximated initially by the homogeneous isotropic 
Friedmann- Robertson- Walker form in (JHJ) with (3i ~ and = dx % . If it? > 1 and the 
p-form terms are negligible, Q a — > as the crunch approaches. For w = 1, f2 CT remains small 
but finite. If the isotropic case is unstable, then adding anisotropy cannot restore stability; 



just as in section Hi Bl the isotropic scale factors are the most stable. 

It can be shown that the p-form terms involving the spatial gradients of F grow slower 
than the leading homogeneous time-derivative terms, another example of the ultralocal 
behavior discussed previously. Hence, we neglect all spatial derivatives of the field strength. 

The components of F with purely spatial indices, Fi v ..i p+1 are called magnetic and the 
components with one time index, F 0ll '" lp , are called electric, in analogy with the Maxwell 
action. We will use the labels E and B to indicate their respective contributions. F has a 
vanishing exterior derivative dF = 0. In coordinate notation, neglecting the spatial deriva- 
tives of F, this corresponds to 

d [0 F iv .. ip+l] = 0, (30) 

where the brackets [■ • ■ ] indicate antisymmetrization. Thus, the magnetic components are 
constant, 

F il ... ip+1 = (constant) (31) 

The equation of motion for F is 

V^(e A0 F W2 '"^+ 1 ) = ^(e^i?*"^!) + p/^gA^/ta-Mp+i = o. (32) 

Only one set of Christoffel symbols appears due to the antisymmetry of F. Since = 
Jj: log y/^g and r M M j = 0, we can integrate to find, 



F oh-i p = ^ = x ( constant ) (33) 
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The p-form part of the stress-energy tensor is 

\<j> 

T,u = jj^({p+l)F im ...^F v ^"^ - \g^F 2 ). (34) 

Decomposing ()34|) into electric and magnetic components, and including factors of the met- 
ric, we compute the energy density for the p-forms p p = — T °, which is, 

\4> 

n — (P±] L pOii-i p pOji-jp n „. . ilp . p. . n hji . . . n ip+ij P +i\ ( oz\ 

Pp — _^ j\j i 2 r r y^n y^ P j P + 2 r ^- % p+^ r n-]p+iy y ) y 00 ) 

where the positive constants a 2 E and a 2 B represent the magnitude of the electric and magnetic 
energy, respectively. We can now define a new set of fractional energy densities, 

n </l = p- 1 ($ 2 /2 + V(<j>)), (37a) 

n z = P~ 1 ^Tp) a Ei ( 37b ) 
n B = p- 1 - mr) o? B , (37c) 

e -A0 gA0 

p = ±0 2 + V(4>) + a% + ^ry«|. (37d) 

where is the energy density in the scalar field and Qe arid Qm are the energy densities in 
electric and magnetic modes. We are assuming that the anisotropy is negligible, so Q a ~ 0. 
The solution for a ^-dominated universe with equation of state w is, 



<p = q\n\t\, q=J—— (38) 

and a = \t / t \ 2 ' 3 ( 1+w > . Substituting in (JH§J), two terms in p p may be written as 

p p = a 2 E \t\ PE + a 2 B \t\ PB , (39) 
where pe and pb are called the electric and magnetic exponents, respectively. They are, 



3 1 + w) \/ 3 1 + w) 



3(1 + w) V 3(1 + iu) 

Note that these expressions are invariant under a duality transformation, which takes p 
2 — p, interchanges the electric and magnetic modes, and takes <j) — ► —0. 
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In the Friedmann equation, the scalar field energy density scales as t~ 2 . Consequently, 
— > 1 and Qe,b -> as the universe contracts if both pe and pb are both greater than 
—2. In this case, the p-form contribution is negligible and A is non-chaotic. Alternatively, 
if either p B or p E is less than —2, the respective p-form terms become large and alter the 
dynamics. 

For w = 1, the non-chaotic values of A are 



5/3 < A < p = 



~^2j?><\<^2j?> p=l (42) 



0<A< v / 8/3 p = 2 

Increasing w causes the interval of non-chaotic couplings to grow, as shown in Fig. |21 In 
particular, for any p and A, there exists a critical value w crit (X,p) such that, for w > w crit (X,p) 
the p-form terms remains negligible. For any set of p-forms and couplings there exists a w CT it, 
the maximum of 1 (the critical equation of state for curvature) and the w CI i t (X,p) for each 
p and A. Then the contraction is non-chaotic if w > w CT n. 

The behavior can be understood in terms of an effective equation of state for the action 
(J2BJ), using the conservation equation 

p=-3~(l + w eS )p, (43) 
a 

where p is given by (|37dj) . Using (J37j) . the equation of motion for <fi and the Friedmann 
equation, we find 

Wea = + + ^n B , (44) 

where 

wa = -. (45) 

2 /2 + V(0) 

is the equation of state for the decoupled scalar field and Q^ + Qe + Qb = 1. The expression 
(144)1 is exact, valid for all values of the assuming the background is homogeneous, flat 
and isotropic. For the electric and magnetic contributions, we can introduce We = ^y 2 and 
wb = respectively. The w e g is just the f2- weighted average of w^, we and wb- 

All the A dependence of w e fr is contained in the time evolution of the f^; w^, we and wb 
do not depend on A. Both we and wb are always less than or equal to unity, and at least 
one is strictly less. 
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FIG. 2: The four dimensional electric and magnetic couplings A as a function of the critical equation 
of state for p = 0, 1, 2. The upper and lower three curves represent the critical electric and magnetic 
exponents, respectively. A form with given p and A is stable in a universe with equation of state 
w if the point (w, A) lies between the two curves for the given p. 

If the p-form coupling A is non-chaotic, the behavior is simple. The quantities Qe and 
Qb rapidly approach zero as approaches one, and the universe is dominated by the scalar 
field, with the equation of state w<$,. This is the non-chaotic case, discussed in section ITTT1 

Alternatively, if the p-form coupling is supercritical, Qe and Qb grow. The averaging 
of the E and B component ensures w e s < w^. If — 1 then iy eff < 1. In this case, 
the anisotropy grows and chaotic oscillations occur. It is not known if this happens in the 

> 1 case. If in addition, the p-form coupling is critical (so w cr \t = 1), it turns out that 
the model is equivalent to an infinite-dimensional hyperbolic Toda system. There are an 
infinite number of jumps from one Kasner-like solution to the next, but the system may be 
formally integrable 0,0]- It is not clear what the physical ramifications of this behavior are. 
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IV. TIME- VARYING EQUATION OF STATE AND CHAOS 

In a realistic cosmological model, the equation of state will not be constant, but will 
depend on the scale factor and approach some limiting value w — > w as a — > 0. If w ^ w cr i t , 
none of the above analysis changes substantially. The model is supercritical if w < w crit 
or non-chaotic if w > w cvlt . The critical case, w = w cr i t , is more subtle, and the time 
dependence of w can be significant. In this section, we assume w CTit = 1, as this is the most 
important case, and analyze what happens when — > 1 at the crunch. We can expand 
as 



where 7 is a small function of the scale factor such that 7 — > as a — > 0. 

If there is no p-form with critical coupling, then using (JJJ) and (j!5|) . it can be shown that 
if 7(a) log a approaches a constant as a —>■ 0, then the behavior is essentially the same as 
the w = 1 case, i.e. non-chaotic. The radius of the Kasner circle in figure ^ shrinks, if 
w — ► 1 + , or expands, if w — > 1~. If 7 — > so slowly that 7(a) log a diverges as a — > 0, then 
the anisotropy is eliminated if 7 approaches zero from above or the chaos is restored if 7 
approaches from below. 

Alternatively, if the model has a p-form with critical coupling, the Kasner contraction 
will be stable if the p-form contribution to the equations of motion remain subdominant, or, 
equivalently, if the ratio of the p-form terms to the other terms vanishes in the a — > limit. 
For magnetic modes with critical coupling A cr i t < 0, we find: 



where C is a positive constant and by ~ we mean up to terms finite in the a — > limit. The 
behavior is identical if the electric modes have critical coupling. If 7 — > very slowly, for 
example 



so that the integral diverges as a — > 0, then the ratio goes to zero and VLm becomes negligible 
in the a — > limit. This ensures that the term is small, and never grows to influence the 
dynamics. 

Let us investigate what conditions on the potential will give us a 7 of this form. If we 




(46) 




(47) 




(48) 
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combine the Friedmann equation and equation of motion for 0, we obtain 



where , denotes a derivative by and 

= (50) 
a log a 

The equation of state can be expressed in terms of if;, 

1^ = 1 + 7 = 2^-1. (51) 

We can obtain — > 1 + as a — > for any negative potential which is bounded (for large 
negative values of 0) by — Ce~^^, where C is a positive constant (see (|29|)). The kinetic 
energy increases more rapidly than the potential energy in these cases, and so approaches 
unity at the crunch. In particular, the potential need not be bounded below. In general, 
any potential which can be expressed in the form 

V(<t>) = 2W'((f)) 2 - 3W((j)) 2 (52) 

32] ]. Hertog et al (^J have shown that the potential 

-V e-* (53) 



satisfies positive energy 



where Vo and c are positive constants, can be expressed in this form provided c < a/6, and 
so satisfies positive energy. For c > y/E, solutions exist with total ADM energy that is 
unbounded below. 

For the potential ([53")). V^/V = c. In the case c < \/E, we find 

7 oc a y , (54) 

where y is a positive constant. Consequently, 7 log a — > as a — > and the p-form with 
critical coupling is not suppressed. However, when c = y/E, the solution to the equations 
of motion show that 7 log a approaches a constant, so the p-form can be suppressed when 
positive energy is violated. 
The potential 

V{<j>) = -Voe-^l", (55) 
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(or more generally, an exponential times any finite order polynomial) satisfies positive energy 
(i.e. can be expressed in the form (|52j)) for n < —1. Solving the equation of motion ()49|) 
for large 0, we find that for n < — 1 the p-form with critical coupling is not suppressed. 
Surprisingly, for n > — 1 the ratio (J47|) goes to zero, and the solution is stable. For the 
broad class of potentials (|53|) and (|55|) . the parameters for which they satisfy positive energy 
turn out to be exactly those which do not suppress the p-form. It is an open question 
whether any potential can be constructed which will suppress the w cr i t = 1 p-form and 
satisfy positive energy. 

V. EXTRA DIMENSIONS, ORBIFOLDS AND CHAOS 

In models in which gravity is fundamentally higher dimensional, the detailed global struc- 
ture of the extra dimensions can suppress or enhance chaos in the four dimensional theory. 
We consider two simple compactifications of five dimensional gravity, on S 1 and 5 ,1 /Z 2 . In 
the first, the chaotic nature of pure five dimensional gravity descends to the four dimensional 
theory. In the second, the chaotic behavior is suppressed. Models of quantum gravity also 
generally include additional matter fields in the extra dimensions. As an example, we dis- 
cuss the compactification of heterotic M-theory to four dimensions and find that its behavior 
during gravitational contraction is on the borderline between smooth and chaotic. 

Consider a five-dimensional, flat universe without matter fields. We know from the study 
of general Kasner universes that it will exhibit chaotic behavior. Now compactify one 
dimension on S 1 . We know that the four-dimensional effective theory describes Einstein 
gravity coupled to a free scalar field. The scalar field describes the volume of the S 1 - it is a 
simple example of a moduli field. As all of our preceding arguments regarding gravitational 
contraction are local in nature, we expect that the resulting system should be chaotic as 
well. However, a free scalar field has equation of state w — 1. According to our analysis in 
section m one might think that the behavior should be non-chaotic. What has happened to 
the chaos? 

The resolution lies in the fact that we have neglected many of the degrees of freedom of 
the higher-dimensional theory. A general five dimensional metric Gmn can be written, 




(56) 
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If we neglect the dependence of the metric on the fifth dimension, we may integrate out 
this dimension and perform a conformal transformation to canonically normalize the four 
dimensional Ricci scalar. The coefficient q is then chosen to canonically normalize the scalar 
field kinetic energy in the resulting action 

S = J d A x^{\R - \{dct>) 2 - |e^F 2 ), (57) 

which describes a vector field coupled to a free scalar and to gravity. The coupling A = 
is outside the stable range for a 1-form in four dimensions. Therefore, the four dimensional 
theory is chaotic, as we would have guessed, but we have to include the interactions with 
p-forms to see that this is so. 

Next, instead of compactifying the fifth dimension on S l , let us compactify on the orbifold 
S l /7i2- If the coordinate x 5 on S 1 runs from — 7r to it, this orbifold can be realized as S 1 
together with the reflection x 5 — > —x 5 . This takes — > —G^, or equivalently — > — A^. 
Thus, the Kaluza-Klein zero-mode vector field A^ is absent in the effective action (J57|) . 
The absence of this vector field in the effective action thus implies that the four dimensional 
theory is no longer chaotic. 

While orbifolding suppresses some gauge fields and p-forms that would cause chaotic 
behavior, in some models there are additional p-forms in the bulk. These p-forms, after 
dimensional reduction, may themselves lead to chaotic behavior. An illustrative example is 
heterotic M-theory, which includes a three-form field. 

The low-energy four dimensional effective action has been evaluated perturbatively by 
Lukas et al js^J. To zeroth order in the eleven dimensional gravitational coupling k, it is 

S (0) = ^JV_ J dx 4^( R _ {da) 2 _ {dc) 2 _ e -^T 3c{dx) 2 _ e -V8a ((9a) 2) ? (5g) 

where we have rescaled the fields in Lukas' action so the kinetic energies are canonically 
normalized. The scalar field c is the radion, which governs the brane separation. The 
Calabi-Yau volume modulus a and scalar field a (which comes from the eleven dimensional 
3-form) do not couple to c, and so can be ignored. However, the 3-form modulus \ couples 
to c and the exponent is critical A = —a/8/3. Hence, the theory does not lead to stable 
Kasner contraction. Including the first order (k 2//3 ) correction to the action does not change 
the result. 

As this theory is critical, it is quite conceivable that higher order corrections will lead 
to a different behavior during cosmological contraction. There are a number of kinds of 
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corrections to (|58j) that could push the theory away from criticality and render it either 
chaotic or non-chaotic; but, it is not yet known which behavior occurs. 

VI. CONCLUSIONS 

The new results presented in this paper build on over three decades of preceding research 
on the behavior of cosmological models contracting to a big crunch. The classic work focused 
on cases where the equation of state of the dominant energy component is w < 1 and w is 
constant. The essential results in this case are: 

• For a perfect fluid with w < 1, the contraction is smooth and anisotropic in the absence 
of curvature and chaotic mixmaster if there is non-zero curvature. 

• For a perfect fluid with w = 1, the contraction is smooth and anisotropic in the absence 
of curvature. With curvature, the contraction is anisotropic also, although, depending 
on the initial anisotropy, the contraction may undergo a finite number of jumps from 
one Kasner-like behavior to another. 

• For a free scalar field coupled to p-forms with coupling e A( ^, the contraction is chaotic 
mixmaster if the coupling A is outside a finite interval of non-chaotic A. The mixmaster 
case is non-integrable and the critical case may be integrable. 

In this paper, we have extended this work to include cases where w > 1, a situation that 
arises naturally in some recent models with a big crunch/big bang transition, such as the 
cyclic and ekpyrotic models. We have added the following results: 

• For perfect fluid with w > 1, the contraction is smooth and converges to isotropic at 
the crunch. The Einstein equations converge to ultralocal, homogeneous and isotropic 
Friedmann equations. 

• For a scalar field coupled to p-forms, there exists a w CT i t such that the contraction is 
smooth and isotropic for w > u> cr it- 

• If w is time-varying and approaches one from above sufficiently slowly the contraction 
is smooth and non- chaotic, even in the presence of a p-form with critical equation of 
state w CTit = 1. 
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• In models with an extra dimension, compactification generically produces a scalar field 
and p-forms. Z2 orbifolding forces some p-forms to zero and, thereby, suppresses their 
contributions to chaos. 

In this paper we have studied how chaotic mixmaster behavior may be suppressed in models 
involving a big crunch/big bang transition. In particular, the ekpyrotic and cyclic models 
already include some of the required ingredients including a scalar field with w > 1 and 
Z 2 orbifolding. We did not present a complete non-chaotic string-motivated model, but the 
considerations reported here will we hope be helpful in that regard. 
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